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Abstract. The configuration model is one of the most successful models for generating uncor- 
rected random networks. We analyze its behavior when the expected degree sequence follows a 
power law with exponent smaller than two. In this situation, the resulting network can be viewed as 
a weighted network with non trivial correlations between strength and degree. Our results are tested 
against large scale numerical simulations, finding excellent agreement. 



INTRODUCTION 

Complex networked systems represent better than any other the idea of complexity. 
When looking at their large scale topological properties, real networks are far more 
complex than classical random graphs lUlU, showing emerging properties not obvious at 
the level of their elementary constituents -the small-world effect, scale-free connectivity, 
clustering, degree correlations, etc. When looking at the dynamical processes that take 
place on top of them, these large scale topological properties have striking consequences 
on the behavior of the system -absence of epidemic threshold, resilience to damage, 
etc. The understanding that many real world systems of interacting elements can be 
mapped into graphs or networks has led to a surge of interest in the field of complex 
networks and to the development of a theoretical frameworks able to properly analyze 
them JUmBI]. Under this approach, the elements of the system are mapped into vertices 
whereas interactions among these elements are represented as edges, or links, among 
vertices of the network. 

In an attempt to bring nearer theory and reality, many researchers working on the 
new rapidly evolving science of complex networks have recently shifted focus from 
unweighted graphs to weighted networks U i i, M & 

Commonly, interactions 

between elements in network-representable complex real systems -may they be com- 
munication systems, such as the Internet, or transportation infrastructures, social com- 
munities, biological or biochemical systems- are not of the same magnitude. It seems 
natural that the first more simple representations, where edges between pairs of vertices 
are quantified just as present or absent, give way to more complex ones, where edges are 
no longer in binary states but may stand for connections of different strength. 

In an unweighted network, all the topological properties can be expressed as a func- 
tion of the adjacency matrix, Ay, whose elements take the value 1 when vertices i and j 
are connected by an edge and otherwise. Using this matrix, the degree of a vertex -the 



number of neighbors or connections it has- becomes 

The distribution of P(k), is called the degree distribution of the network, and it 
arises as the most fundamental topological characteristic. Surprisingly, in a vast majority 
of cases, real networks show degree distributions following power laws of the form 
P(k) ~ kr* 1 for k ^> 1 and 2 < 7 < 3. This implies that the second moment of the degree 
distribution, (k 2 ), diverges in the thermodynamic limit, which causes the loss of any 
characteristic degree scale. For this reason, these class of networks are called scale- 
free (SF). This feature of the degree distribution is, eventually, the responsible for the 
surprising behavior of dynamical processes that run on top of these networks. 

For weighted networks, the adjacency matrix is no longer the fundamental quantity 
ruling the properties of the network. Instead, we have to consider the matrix coy, which 
measures the weight between the pair of vertices i and j, that is, the magnitude of the 
connection between i and j. The analogous quantity to the vertex degree is now the 
vertex strength, defined as 

Si = Y i ( °ii A H- (2) 

j 

This quantity measures the total strength of vertex i as the sum of the weights of all its 
connections. If the weights of a given vertex are not correlated with the vertex degree, 
the average strength of a vertex of degree k t is simply given by 

s(ki) = (a)k u (3) 

where (co) is the average weight of the edges of the network. In this situation, strength 
and degree are proportional and weights do not incorporate more information to the 
network than that already present in the adjacency matrix. Real networks show, however, 
a very different scaling, with a non trivial relation between strength and degree of the 
form 

(4) 

with j8 1 (typically /3 > 1). This anomalous scaling reveals the presence of correla- 
tions between strength and degree and, thus, the need to model network formation in 
an weighted basis, where the evolution of the network is linked to the evolution of the 
weights assigned to connections. In this paper, we present a model of weighted net- 
work based on the configuration model. We shall show that, when the expected degree 
sequence follows a power law with exponent smaller than 2, the resulting network is 
weighted and shows a non trivial scaling between strength and degree. 



THE CONFIGURATION MODEL 



The configuration model was first introduced as an algorithm to g enerate uncorrected 
random networks with a given degree distribution, P{k) lfl2l Il3l IbH The model 



operates in two steps: 

• We start assigning to each vertex i, out of a set of N, a number of "stubs", drawn 
from the distribution P(k), under the constraint that the sum is even. 

• Pairs of these stubs are chosen uniformly at random and the corresponding vertices 
are connected by an undirected edge. 

Given the random nature of the edge assignment, this algorithm generates networks with 
the expected degree distribution and no correlations between the degrees of connected 
vertices. The model, as described above, allows the formation of multiple or self connec- 
tions among vertices. Nevertheless, when the expected degree distribution has bounded 
fluctuations, the number of such "pathological" connections is small an can be neglected 
in the thermodynamic limit. In this case, one can add an extra constraint in step two of 
the algorithm avoiding multiple or self connections without modifying the resulting net- 
work. 

The picture changes drastically when the expected degree distribution has unbounded 
fluctuations. This is precisely the case of SF distributions with exponent 7 G (2,3]. In 
this situation, vertices with expected degree satisfying ki > ■sj {k)N cannot avoid to 
form multiple connections. Indeed, as it was proved in [16], the number of multiple 
connections in this case scales as N 3 ~ r lnN. Yet the situation is not so dramatic since the 
overall number of connections is much larger that the number of multiple connections. 
However, one has to be careful because the "hubs" of the network -those in the tail of 
the distribution- are precisely the ones more prone to hold multiple connections and, 
therefore, it could alter the results of dynamical processes evaluated on top of these 
class of networks. There are two strategies one can follow in order to avoid multiple 
connections: 

• One can introduce a constraint in step two prohibiting multiple connections. This 
has the side effect of introducing strong disassortative degree correlations among 
connected vertices. 

• Alternatively, one can introduce a cut-off in the distribution P{k) scaling as y/N. 
By doing so,one recovers uncorrelated networks but with a smaller cut-off than the 
natural one IU7I1 . 

THE CONFIGURATION MODEL BEYOND 7=2 

The configuration model with expected degree distribution with 7 < 2 has not been stud- 
ied so forth. This is an extreme situation in which the number of multiple connections 
cannot be neglected. We can take advantage of this fact to construct a weighted network, 
where the weight between a pair of vertices is the number of multiple connections they 
have lfllL[l8ll . Now, the expected quantity of a given vertex is no longer its degree but 
its strength. Therefore, let us change notation and denote by P(s) ~ s~( 1+T \ < z < 1, 
the distribution of expected strengths, that is, the number of expected connections of 
vertices. 



At the level of strength, vertices are uncorrelated and the average weight of edges 
linking nodes of strengths Si and Sj can be written as 



% = (5) 

where S is the total number of connections -regardless they are multiple or not- in the 
network, that is, 

S= l -£ Si ~N^ (6) 

i 

The average degree of a vertex can be obtained as 

k = £ tiHj®(l - GHj) + £ ®{GHj - 1) , (7) 

j j 

where ©(■) is the Heaviside step function. The first term in this equation is the con- 
tribution of those connections that, on average, are smaller than 1. The second term 
represents connections with an average number larger than 1 and, thus, holding multiple 
connections. In this case, the contribution to the degree is just 1 and not Ctyj. Using the 
expression for the weights, Eq. ©, in the continuum limit we can write 



r2S/sj s . s . roo 

k( Si )=N / p( Sj )±±dsj + N / P(sj)d Sj (8) 

J\ m J2S/si 

that, in the thermodynamic limit goes as 

It is worth noticing that there are finite size effects that depend on the specific form of 
the strength distribution. The model generates a non trivial weighted network with an 
exponent /3 = 1/t larger than 1, as it is observed in real weighted networks H2D - 

The tail of the degree distribution can be obtained from this scaling relation, yielding 
an exponent 7= 2. Thus, even though the expected strength distribution has an exponent 
smaller than two, the resulting network has exponent equal to two and a number of 
connections growing as Y,ih ~ NlnN. 



NUMERICAL SIMULATIONS 

To test the results of our analysis, we have performed numerical simulations of the 
configuration model, generating expected strength distributions of the form 

p «~(d^- (10) 



The value of s c modulates the probability to find strengths of value larger than 1 . We 
choose this particular form because, by an appropriate choice of s c , the finite size effects 
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FIGURE 2. Simulation results for the average degree of vertices of strength s, k(s), as a function of 
s, for T = 0.95, T = 0.85, and T = 0.75. The size of the network is N — 10 5 . In all cases, the strength 
distribution is given by Eq. i ll Oi l with s c = T. Solid lines correspond to the best fit estimates (numbers in 
parenthesis into the legend box) of the theoretical behavior k(s) ~ s T . 
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FIGURE 3. Cumulative degree distributions, P c (k) = Y^k'=kP(^)' generated by de model. The solid line 
is power law of the form k~ l corresponding to an exponent 7=2. 



are minimized. In all the simulations, the size of the network is N = 10 5 , the values of 
x are 1.95, 1.85, and 1.75 and s c = x. We first show, in Fig.d the cumulative strength 
distributions used in the simulations as compared to the theoretical ones, 

(s — xy 

computed using the continuum approximation. 

Fig. shows the scaling relation between k(s) and the strength s. As it can be 
seen, the theoretical prediction s T is well satisfied, with the following estimates for the 
exponent JS -1 : j8 1 = 0.74 ±0.01 for x = 0.75, /3 1 = 0.84 ±0.01 for x = 0.85, and 
/3 1 = 0.95 ±0.01 for x = 0.95. Finally, the resulting cumulative degree distribution is 
plotted in Fig. El showing that this function goes, for large degrees, as k^ 1 independent 
of T, as predicted by our analysis. 



CONCLUSIONS 

In summary, we have analyzed the behavior of the classical configuration model in the 
case of expected strength distributions following a power law form of exponent smaller 
than two. We have shown that, in this case, the resulting network is weighted, where 
the weight stands for the number of multiple connections among vertices. The model 
presents a non trivial scaling relation between strength and degree and a degree distri- 
bution with exponent 7=2, independent of the exponent of the strength distribution, x. 



These results highlight the subtleties that may arise when dealing with SF networks even 
in the most simplified models. 
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